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BACKGROUND OF THE INVENTION 

I. Cross Reference to Related Materials 

5 Previously assigned patents contain information relating to the determination of certain physical parame- 

ters with respect to fuels of interest in the present application. These include: U.S. Patent 4 944 035 entitled 
MEASUREMENT OF THERMAL CONDUCTIVITY IN SPECIFIC HEAT which issued on July 24, 1990 to Roger 
L Aagard, Ulrich Bonne, the inventor in the present application, and Robert J. Matthys; U.S. Patent 4 956 793 
which issued September 11, 1990 to Ulrich Bonne, the inventor in the present application, and Steven D. 

10 James; U.S. Patent 4 961 348 issued October 9, 1 990 to Ulrich Bonne, the inventor in the present application, 
entitled FLOW METER FLUID COMPOSITION CORRECTION; Serial No. 07/285,890 filed December 1 6, 1 988 
entitled LAMINARIZED FLOW METER; Application Serial No. 07/789,411 filed November 1, 1991, which is a 
continuation of application Serial No. 07/429,138 filed October 30, 1989 to Ulrich Bonne, the inventor in the 
present application. To the extent necessary for the complete description of any aspect of the present appli- 

15 cation, material from the above may be deemed incorporated by reference into this application. 

II. Field of the Invention 

The present invention is directed generally to combustion control and, more particularly, to the improved 
20 determination of fuel characteristics including heating value and oxygen demand and Wobbe Index or Number 
which can be utilized for more accurate fuel metering and combustion control (i.e. efficiency and firing rate) in 
terms of accurately determining the amount of oxygen needed for complete combustion for a given type and 
amount of fuel which is valid for a wide range of fuel compositions and more accurate inasmuch as it is less 
sensitive to errors in the measured parameters from which the values are derived. 



25 



III. Description of the Related Art 



The determination of the oxygen demand and heating value, energy content or heat of combustion of a 
fuel represents an important determination with respect to both the operation of a combustion system with 

30 respect to proper combustion of the fuel and in determining the quantitative commercial value of the fuel with 
respect to the supplier. In this regard, various methods for measuring the heat available from fuel gases have 
been implemented and used for many purposes. As evidenced by the above-referenced patents, the use of 
spaced microscopic heating and sensing elements on semiconductor chips to measure certain physical para- 
meters in gaseous media is known. These "microbridge" system are extremely fast reacting, very stable and 

35 very sensitive with respect to measuring thermally induced changes in electrical resistance when placed in a 
fluid of interest. The use of such devices to accurately determine thermal conductivity, k, and specific heat 
are shown in Aagard, et a!., above, and Bonne, et al. 793, also cited above, also enables the accurate deter- 
mination of specific gravity of a species of interest 

The related application, Serial No. 07,789,411 filed November 1, 1991, discloses the use of k and Cp for a 

40 determination of heating value or heat of combustion, H, which is, in turn, used for the determination and control 
of the delivered energy and the oxygen demand to control the efficiency or fuel-to-air ratio in a metering and/or 
combustion system, respectively. All determinations of Y = H or D or Wo for a given type and amount of fuel 
or the oxygen demand of the fuel, D, were correlated to k and Cp. The relationship utilized in the above cross- 
referenced application Serial No. 07/789,411 was written in terms of a direct relation of thermal conductivity, 

45 k, and specific heat, Cp, as measured for flowing fuel and evaluated according to the following polynomial re- 
lationship: 

H = A^'Mx) • A 2 f 2 "2(x) ■ Aafa^x) (1) 

where 

A 1f A2, A 3 = constant coefficients 
so n1, n2, n3 = exponents 

f,(x) = k (thermal conductivity at a first temperature) 

hi*) - k (thermal conductivity at a second temperature) 

Wx) = Cp (specific heat determined at one temperature, normally the temperature of f t (x) or f 2 (x)). 

Results obtained utilizing the relationship of equation 1 for a selection of over 60 natural gases was used 
55 to obtain empirical constants and exponents, and these, in turn, were applied to later determinations of different 
fuels. The results obtained utilizing equation 1 with the empirically derived constants and exponents have 
proved to be generally quite satisfactory calculated (error generally was < 0.5%) for natural gases characterized 
by concentrations of carbon dioxide (COJ and nitrogen (NJ not exceeding two percent (2%) by volume. 
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When the same equation was recently applied to an expanded range of gases notably including gases con- 
taining nitrogen levels as high as over 20% and containing much higher concentrations of ethane, propane, 
C0 2 and hydrogen, however, errors in the determination of heating value increased from less than 0.5% to 
over 1 2%. This is illustrated for 78 natural and 22 test gases in Figure 3. Even after extensive work in readjusting 
5 the empirical coefficients and exponents of the algorithms of the equation 1, the maximum error has been > 
found to be still over 6%. An error this large still is unsatisfactory for most uses of the higher heating value, 
for the purposes of substantiating amounts charged customers for the heating value of fuel supplied, or fuel- 
to-air ratio combustion control based on oxygen demand measurement In this regard, any improvement which 
substantially reduces the error and increases the reliabOity of determinations of the type involved without add- 
to ing a great deal of cost to the measurement apparatus would be highly desirable. 

The invention presents a highly accurate on-line system for the determination of the higher heating value 
H by monitoring the heat content of gaseous hydrocarbon fuels on a volumetric basis. In addition, however, 
knowledge of another very important dimension is necessary to account for changes in the specific gravity of 
the fuel so that a constant heat input to a pressure fueled burner system can be maintained. Thus, for example, 
15 while utility companies are very much concerned that the customer be charged precisely for the amount of 
heat value sold on a volumetric basis, those utilizing the gaseous fuel are also concerned that the Wobbe Num- 
ber of Index, defined below, be considered so the burner energy input remains constant, whatever the relative 
Btu content be on a volumetric basis. This is especially important in devices that use mass flow devices such 
as orifice meters to measure fuel flow to the burners. Variations in burner efficiency can be very costly espe- 
20 dally in industrial burner applications for which we need to know oxygen demand of the fuel in order to control 
the fuel/air ratio. 

Accordingly, changes in the density or specific gravity of the fuel should be considered in addition to the 
heating value of the fuel on a volumetric flow basis in order to assure constant heat input to the burners. Such 
may be accomplished by determining and monitoring the Wobbe Number or Index of the fuel, Wo, which is 
25 defined as follows: 

Wo = Hp - o- 6 (2) 

where: 

H = the higher heating value of the fuel in Btu'stft 3 joules/meter 3 , etc. and 
p = the relative density or specific gravity of the fuel. 

30 It is apparent from the above that by considering changes in the specific gravity of the fuel gas in addition 
to the heating value on a standard volumetric basis, the Wobbe Index is useful as a measure of the potential 
heat production available for a given burner input Thus, a fuel mixture of constant Wobbe Index will provide 
constant heat input to an orifice-controlled burner system. Because of the great variation in the heating value 
of gaseous fuels, the Wobbe Index has long been used in connection with providing a more constant firing 

35 rate input to burner systems to produce the specified performance and system stability. In the past, however, 
in order to determine the Wobbe Index of a fuel gas, separate calorimetry and density measurements had to 
be made and the results combined. This involved the use of expensive calorimeters and density meters and 
required considerable time to obtain useable results. It was generally not affordable to determine a represen- 
tative real-time Wobbe Index on-line control in an industrial setting. 

40 For diverse applications, the measurement or determination of thermophysical or thermochemicaf related 
properties is very desirable. These include compressibility (Z), viscosity fa), pseudocritical temperature and 
pressure. Other parameters or fuel properties of interest include f lammability limits, critical compression ratio 
(i.e. the maximum compression ratio before incipient knock), flame speed, and yellow tipping of the flame. Of 
course, it would be very desirable if all of these parameters could also be accurately correctable to the meas- 

45 unable or readily determinable quantities k, Cp, dk/dT, dCp/dT. T and P. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, there is provided an improved method for the determination of 
so many thermophysical and thermochemicaJ parameters of fuel gases including but not limited to heating value 
or heat of combustion, H, amount of oxygen needed for complete combustion for a given type and amount of 
fuel, or oxygen demand, D, and additional parameters as the Wobbe Number or Index, Wo, based on new dis- 
coveries with regard to determining the relationship among the above and the thermal conductivity, k, and 
specific heat, Cp. The present invention contemplates the use of a class of relationships among these values 
55 which has been found to be valid and more accurate, for a much wider range of fuel compositions including 
those with large amounts of higher alkanes, hydrogen and non-combustibles. The present invention provides 
this accuracy in a manner which exhibits a lower sensitivity to measurement errors in obtaining k and Cp and 
which is capable of readily presenting yield outputs which can be utilized for more than one application or usage. 

3 
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Unlike previous determinations, the present invention utilizes parameters which are directly measurable are 
first converted to a reference condition of pressure and temperature and may not require the measurement of 
pressure. This has led to a new and more accurate polynomial expression of H, Wo, or D. The new, more ac- 
curate polynomial representation uses a structure leading to more accurate results. 

As indicated above, while earlier expressions of H or D in terms of measured values of k or Cp were accurate 
(i.e. < 0.5% error) with respect to pure fuels or fuels having small amounts, i.e. less than about 2% of non- 
combustibles such as C0 2 or N 2 at higher nitrogen levels, i.e. 20% or more, or higher concentrations of ethane, 
propane, C0 2 , H 2 , or the like, the accuracy of such determinations were subject to a great deal more error (> 
6% even after readjustment of the coefficients and exponents of the algorithm of equation 1). 

According to the present invention, it has been found that the relationship that exists between H, Wo or 
D, and k and c p can be expanded to accurately cover a greater variety of fuel gases, including those containing 
large amounts of alkane hydrocarbons, higher than methane outside the main range of conventional fuels and 
containing higher amounts of other species including C0 2 and H 2 without sacrificing accuracy of determi- 
nations obtained with respect to the other fuel gases. In addition to greatly expanding the variety of accurately 
determinable fuel gases, the present invention also makes it readily possible to measure or determine other 
thermophysical or thermochemical related properties including compressibility (Z), viscosity ft), pseudocriti- 
cal temperature <Tpc) and pseudocritical pressure (PpJ, flammability limits, critical compression ratio, flame 
speed and yellow tipping of the flame. 

According to the present invention, it has been discovered that a generally universal correlation exists be- 
tween Y, a given property of interest, and the measured microbridge properties of thermal conductivity, k, spe- 
cific heat, c p , temperature, T, pressure, P (which can also be independently sensed), the derivatives, dk/dT 
and dCp/dT, and ratios ki/k 2 , Cp^ and c^/Cp^. 

A preferred relationship in accordance with the invention can be expressed as follows: 



Y s = A 0 + A^nlc^l x s Ply s ql T s rlp s Sl + 
+ A 2 k S n2 ^p m2 x s P2y s q2 Ts r2p s s2 * . . . 



= EiAikgHic^mi x s Piy s qi V^si (4) 

where 

dependent variable, Y„ represents any one of higher heating value, H; oxygen demand, D; Wobbe Index, 
Wo; relative density or specific gravity, p; absolute density Pa ; inerts, I; compressibility factor, Z; critical com- 
pression ratio, Rc; viscosity, x\; etc. 

A©, A, . . . A, are constants or coefficients, 

ni, m,, p,, q } , r,, s, are exponents of values from 20 to 20, including zero, with the number of terms, i, ranging 
from 1 to 15, 

kg and represent their values at a predetermined condition, which could be 60°F, 0°C, 15°C, or any 
other arbitrarily chosen temperature and pressure, 
Xg represents dk/dT at the reference condition, 

y s represents dCp/dT or c^ Ti Ic^ T9 at a reference condition, T, * T s (the units of Cp are typically energy 
per unit of mass or mole and per AT, while those of Cp, are typically energy per volume and per AT), 
T s represents temperature at said reference condition, 
P s represents absolute pressure at the reference condition. 
A preferred method of determining kg and Cp, may be represented by the expressions: 
kg = ao + atkfl"* - (k* - k„) (T, - T 8 )/(T 2 - T,) (5) 
Cps = bo + b, V 1 - Me,* - Cpt,) (T, - T,)/(T 2 - T,) (6) 

where 

a, and b, are constants, 
ml and n1 are exponents 

T, and T 2 are temperatures or pressures In a desired range 

T 8 is the standard temperature or pressure. 
This relation has been found to significantly reduce errors in k and Cp as follows: 

1) It uses values of at least two measurements of k and Cp at different temperatures interpolated or ex- 
trapolated to a standard condition increase accuracy. 
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2) Using the ratio of values at different thermodynamic conditions allows one to make use of physical prop- 
erties that are largely or more independent of such thermodynamic condition (e.g. varying the temperature 
at constant pressure for c^n/c^n eliminates the pressure influence of Cpv whereas it would not be elim- 
inated by simply using dCpv/dT). 
5 3) Even more accuracy may be available if the relation is fitted individually to a series of limited ranges of 

parameter values which might include two or more ranges for any of the involved parameters; for example, 
in a simple case of two ranges of fuel gas heating values one might select H > 1050 Btu/ft 3 and H < 1050 
Btu/fP. 

Thus, the Y 8 polynomial can then be optimized for each given range. Depending on the application, any number 
10 of relatively limited ranges can be implemented. 

As more fully described in the above-referenced Bonne, et al. (U.S. 4 956 793) patent, the relative density, 
or specif ic gravity, p, in relation to that of air. has been found to be a function of Cp and k according to an empirical 
polynomial relationship of the form: 



15 
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* (Cpkr) = a 6 + a 7 c pt? n7 + 



+ a 8 k T8 n8 + m m m (?) 



where 

20 a* . . . a m are constants 

kn-kTm are thermal conductivities at subscript temperatures T r T m and 
n r n m are exponents. 

The Wobbe Number can also be derived from the higher heating value, H, and the relative density or spe- 
cific gravity, p, according to the relation given above in equation (2) or Wo = Hp 0 - 5 . 



BRIEF DESCRIPTION OF THE DRAWINGS 



FIGURE 1 is a block diagram of a combustion system; 

FIGURE 2 is a block diagram of the sensing and signal processing system of the invention; 
30 FIGURE 3 is a graphical representation of a comparison between the actual and computed heating value 
of 78 gases (o) and 22 test gases (A-V) using a prior correlating relationship; 

FIGURE 4 is a graphical representation of a comparison between the actual and computed heating value 
of the 78 natural gases (o) and 22 test gases (A-V) of Figure 3 using the correlating relationship of the 
present invention; 

35 FIGURE 5 is a graphical representation of the error in the determination of the compressibility factor for 
a range of gases using the relationship of equation 14 over a range of temperatures and pressures; 
FIGURE 6 is a graphical representation of the error in the determination of the actual viscosity of the 78 
natural gases based on the relationship of the invention; and 

FIGURE 7 is a graphical representation of actual excess oxygen control performance at three set points 
40 for three different fuel gases over a wide range of burner firing rates based on the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 shows a block diagram of a typical combustion system generally at 10. Heating system 10 includes 
45 combustion chamber 12, fuel valves 14, air blower 16 and combustion controller 18. Fuel enters combustion 
chamber 12 through fuel conduit 20 where it is combined with air blown from air blower 16. The fuel and air 
mixture is ignited in combustion chamber 12 and resulting flue gases exit combustion chamber 1 2 through flue 

22. 

Combustion controller 1 8 controls the fuel-to-air mixture in combustion chamber 1 2 by opening and closing 
so fuel valves 14 and by opening and dosing air dampers in air conduit 17. Combustion controller 18 controls the 
fuel-to- air mixture based on control inputs entered by a heating system operator as well as sensor inputs re- 
ceived from sensors 24 and 26 in fuel conduit 20, and sensor 28 in air conduit 17. 

Sensors 24 and 28 are typically microbridge or microanemometer sensors which communicate with flowing 
fuel in fuel conduit 20 and flowing air in air conduit 17. This type of sensor is described in more detail in the 
55 above-referenced application Serial No. 285,890. 

Sensors 24 and 28 are directly exposed to the stream of fluid flowing past them in conduits 20 and 17, 
respectively. Sensors 24 and 28 are used to directly measure dynamic fluid flow characteristics of the respec- 
tive fluids. Sensor 26, while in contact with the fuel gas, is recessed in a dead-ended cavity and not exposed 
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to direct flow. 

Figure 2 shows a block diagram of the sensing and signal processing system of the invention which may 
be associated with the combustion system of Figure 1. It includes the dynamic or exposed microbridge flow 
sensors 24 and 28 together with static microbridge 26. In addition, an optional pressure sensor 30 and a tenv 

5 perature sensor 32 are provided. The output of these devices is fed as an input to an analog to digital (A/D) 
converter 34 which provides input to a data processing device such as a data microprocessor 36. Display and 
output devices which may take any form which would occur to those skilled in the art is shown at 38. A power 
supply for the system is depicted by 40. 

Microbridge sensor 26 enables other parameters of the fuel to be measured simultaneously with the dy- 

10 namic flow. Sensor 26 can be used for the direct measurement of thermal conductivity, k, and specific heat, 
Cp, in accordance with a technique which allows the accurate determination of both properties. That technique 
contemplates generating an energy or temperature pulse in one or more heater elements disposed in and close- 
ly coupled to the fluid medium in conduit 20. Characteristic values of k and Cp of the fluid in conduit 20 then 
cause corresponding changes in the time variable temperature response of the heater to the temperature 

15 pulse. Under relatively static fluid flow conditions this, in turn, induces corresponding changes in the time va- 
riable response of more temperature responsive sensors coupled to the heater principally via the fluid medium 
in conduit 20. 

The thermal pulse need be only of sufficient duration that the heater achieve a substantially steady-state 
temperature for a short time. Such a system of determining thermal conductivity, k, and specific heat, Cp, is 
20 described in greater detail in above- referenced patents 4 956 793 and 4 944 035. 

It has also been found that once the specific heat and thermal conductivity of the fluid have been deter- 
mined, they can be used to determine the density or specif ic gravity of the fluid. As described above, this tech- 
nique is more specifically illustrated and described in U.S. Patent 4 956 793. Of course, these parameters can 
be determined by other means if such are desirable in other applications. 
25 Once k and Cp or Cpv of a gas or fluid are known, flow correction factors in the form of simple, constant 
flow-independent factors for the fuel can be calculated and used. The flow correction factors have been de- 
veloped to compensate mass or volumetric flow measurements for changes in fluid temperature, pressure, 
and/or composition. In other words, once k and Cp of a fluid or gas or fuel is known, its true volumetric, mass 
and energy flows can be determined via the corrections: 
30 G * = G(k/ko) m1 (Cp/Cpo) 1 " 2 (Pr/Pro)™ 3 flTTo)"* (7) 

V * = V(k/ko) n1 (Cp/Cpo)"* (Pr/Pro)" 3 (T/To) 04 (8) 
M * = M(k/ko)P 1 (Cp/Cpo)p2 (Pr/Pro)** (T/To)* (9) 
E * = E(k/^Mc^c^ (Pr/Pro)^ (T/To)^ (10) 
where the subscript V refers to the calibration condition such as methane at T 0 and Do and the m, n, 
35 p, and q are optimized but constant exponents; P f in the Prandle Number and G* equals the corrected value 
of the sensor signal G, V* equals the corrected value for the standard volumetric flow V, M* equals the cor- 
rected value for the mass flow, and E* equals the corrected value for the energy flow, E. This technique of 
applying correction factors to the sensor signal, the mass flow, the volumetric flow and the energy flow are 
explained in greater detail in U.S. Patent4 941 348 and co-pending application Serial No. 07/789,411 and which 
40 to any extent necessary is deemed incorporated by reference herein. 

In the present invention, after thermal conductivity, k, and specific heat, Cp, have been determined for the 
fuel flowing through conduit 20, the independent variable, which may represent the heating value, H, of the 
fuel flowing through conduit 20, or other desired parameter, is determined by evaluating the polynomial of the 
form of equations 3 or 4 as follows: 

Y s = A 0 + A 1 k s nl Cps "'l x s Ply s ql T S "P S »1 + (3) 
+ A 2 k s n2 Cp m2 x s P2y s q2 Ts r2p s s2 + . . . 

50 

=S i A i k s ni Cps »i x s Piy s qi T s "P s si {4) 
where 

55 dependent variable, Y 8 , represents any one of higher heating value, H; oxygen demand, D; Wobbe Index, 

Wo; relative density or specific gravity, p; absolute density p a ; inerts, I; compressibility factor, Z; critical com- 
pression ratio, Rc; viscosity, r\; etc. 

Aq, A t . . . A| are constants or coefficients, 
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n». m, f p,, q h r,, s, are exponents of values from - 20 to 20, including zero, with the number of terms, i, 
ranging from 1 to 15, 

ks and Cpg represent their values at a predetermined condition, which could be 60°F, 0°C, 15°C, or any 
other arbitrarily chosen temperature and pressure, 
5 x« represents dk/dT at the reference condition, 

y 8 represents dCp/dT or c^ u fc^ T9 at a reference condition. T, * T 8 (the units of Cp are typically energy 
per unit of mass or mole and per AT, while those of c^ are typically energy per volume and per AT), 
T 8 represents temperature at said reference condition, 
P 8 represents absolute pressure at the reference condition. 
10 in which the needed values of and c^ are determined as in 

K, = a* + a,krf»i - ^(ka - k^t, - y/(t 2 - t,) (5) 
Cps = b 0 + D1 v 1 - Hcpa - Cp,,) (t, - y/(t 2 - t,) (6) 
where a and are constants 
ml and n1 are exponents 
15 t, and t 2 are temperatures or pressures in a desired range 

t„ is the standard temperature or pressure. 
The use of equations (3) and (4) has been found to be very effective in obtaining improved accuracy in 
determining H or D, even outside of the range of normally provided natural fuel gases. A comparison between 
actual and computed heating values of 78 natural and 22 test gases using the relationship of equation (3) or 
20 (4) to determine H is illustrated in Figure 4 in which points A-V represent the 22 test gases and the other points, 
the 78 natural gases tested. In that comparison, the standard computation or aogorithm error, E A , was 2.1 
Btu/f t 3 or 0.21 %. This includes such combinations as ethane or propane and air which have been used as test 
gases or as "peak shaving" gases when the demand for pipeline gas exceeds the supply. In practice, of course, 
the reliability or total error in H or D is not only influenced by the error of the algorithm, E A , but also by its 
25 sensitivity, S, to experimental errors in the input parameters k, Cp, and their temperature derivatives. The im- 
proved algorithm can still result in large total errors, E T , if k and Cp are not measurable with small experimental 
errors, E x , e.g. below 0.1%, because the total error is given by 

E T = SxE x + E A . (11) 

Acomparison between the actual measured and computed heating values of 78 natural and 22 test gases yield- 
30 ed a standard computed error of 1.88 Btu/f t 3 or « 0.1 8%. 

It is well known that hydrocarbon-type fuels combine with oxygen under a constant (hydrocarbon-inde- 
pendent) rate of heat release. The heat released by combustion is about 100 Btu/f t 3 of reacted 0 2 of air at 760 
mmHg and 20°C or (68 °F). This is exactly true for fuel with an atomic hydrogen/carbon ratio of 2.8 and a heating 
value of 21300 Btu/Jb of combustibles and is true to within an error of less than +/- 0.20% for other alkane 
35 hydrocarbons from methane to propane (i.e. CH 4 , C^He and n-Cy-le). 

For example, in order to achieve stoichiometric (zero excess air) combustion, the mixture would be one 
cubic foot of air for each 100 Btu of fuel (e.g. about 0.1 cubic foot of CH 4 ). A more typical mix would be 10% 
to 30% excess air which would require 1.1 to 1.3 cubic feet of air for each 100 Btu of fuel. This would be a 
typical mixture because residential appliances typically operate in the 40-100% excess air range while most 
40 commercial combustion units operate between 10 and 50% excess air. 

Although the present invention has been described with reference to fuels with hydrocarbon constituents, 
the present invention can be used to control the fuel-to-air ratio for other fuels. Each fuel used in combustion 
requires or demands a certain amount of oxygen for complete and efficient combustion (i.e. little or no fuel or 
oxygen remaining after combustion). The amount of oxygen required by each fuel is the oxygen demand value 
45 Dfor that fuel. D, then, is defined as units of moles of 0 2 needed by each mole of fuel for complete combustion. 
For example, the 0 2 demand for CH 4 , C^, C 3 Hb. CO, H 2 and N 2 is Of = 2, 3.5, 5.0, 0.5, 0.5 and 0 respectively. 

Air is used to supply the oxygen demand of the fuel during combustion. In other words, fuel is an oxygen 
consumer and air is an oxygen supplier or donator during combustion. The 0 2 donation Do, is defined as the 
number of moles of provided by each mole of air. The single largest factor which influences D 0 is the humidity 
so content of the air. Absolutely dry air has a value of D 0 = 0.209, while normal room temperature air with 30% 
relative humidity (or 1 volume % or 0.01 mole fraction of H^) has a value of D 0 = 0.207. 

Therefore, one type of combustion, efficiency or fuel/air ratio control can be accomplished by: 1) correlating 
the sensed k and Cp of the fuel to its oxygen demand value D, 2) correlating a second set of sensed k and Cp 
of the air to its oxygen supply, 3) multiplying Dq/D by the desired excess air factor, e.g. A = 1.3 (30% excess 
55 air), and 4) comparing the (set point) AD</D to the measured ratio of standard or actual c but equal T, P of fuel 
and air) volume flow and 5) adjusting and/or to match ADo/D. 

Excess air control performance of a burner in terms of percent excess oxygen based on the relationship 
of the invention is shown in Figure 7 for the three different fuel gases at each of three different excess oxygen 
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set points. The illustrated results represent control parameters that were not corrected for relative humidity of 
air. The burners were operated over a relatively wide range of firing rate ratio, i.e. from 160 MBtu/h to 50 MBtu/h 
(turndown ratio >3:1 ) and maintained, for the most part, on excess oxygen correlation within 1 % or less of the 
set point As seen in Figure 1, the firing rate responds to the demand for heat or steam; the sensors 26 and 
24 determine the or air flow demand to satisfy the set excess air or 02, and the air flow sensor 28 verifies 
that the blower 16 is providing that demand. 

Heating value, H, or Wobbe Number, Wo, would typically be used in fuel for management systems either 
to achieve fuel gases within a desired range of Wo or to achieve an accurate method of billing customers in 
view of the variability of the changed make up of the supply of fuel gases which results in variations of heating 
value and Wobbe Index. Having determined the higher heating value, the Wobbe Number or Index Wo can 
be derived from the higher heating value, H, and the density, p, according to the relation given above in equa- 
tion (2) or Wo = H(p) 0 - 5 . 

The determination of p in absolute density units, as Pa , rather than as specific gravity in relation to air is 
also possible. The actual molar volume V M of the fuel of interest at the measurement temperature and pressure 
is determined based on 

v m = VCpv (12) 

This, in turn, is related to that at a standard temperature and pressure (0°C. 1 Atmosphere), V M0 and the actual 
gas density p a is determined from 

Pa = PoVmo/Vm (13) 

where 

Po is a constant equal to the density of the air at standard conditions (0°C and 1 Atmosphere), 
p A is the absolute specif ic gravity, 

V M0 is a constant equal to the molar volume at standard temperature air pressure, i.e. * 22.4 l/m and 
V M is the actual molar volume at the temperature and pressure of measurement. 
With respect to t he measurement of the compressibility factor, Z, values can either be determined directly 
as a function of of k, Cp, x, y, . . . (see equation 4) or by first computing the pseudo-critical pressure and tem- 
perature, Ppc and Tpe, via correlations, as described below. The known Redlich-Kwong equation may then be 
applied to find the compressibility, Z, by a few iterations: 

Z = 1 + A - AB/(A + Z) + A2B/(AZ + Z2) (14) 

where 

A = 0.0867 P/T r , 
B = 4.934/V-* 
and 

p r = P/Ppc, T r =T/Tpc (reduced pressure and temperature) 
Figure 5 illustrates the correlation and error of this microbridge- based approach. The simulated microbridge 
sensor values of were obtained via equation (14) by using k and Cp values computed for 78 natural gases 
from their known composition which would simulate and represent the result of a microbridge sensor-based 
system (because of the microbridge-measured values of k and Cp); -correct 0 values of Z were obtained by com- 
puting Tpc and via 

Tpo = Wo Zpc = Wc, = ZpqV c (15) 

and 

Ppc = ZpcRoVVp, (16) 
which is also known as the combination rule of Prausnitz and Gunn; the difference 

AZ = ZmP rt mJ^-2P r J r (17) 
is plotted in Figure 5 for the 78 gases, each taken at a matrix of three temperatures and three pressures. As 
shown, the correlation errors at high pressure are largest, but still within about ±0.5%. The total error would 
also be influenced by any experimental errors associated with the k and Cp measurements. 

Viscosity, like compressibility, can also be related to k and Cp by means of equations (3) and (4). This is 
illustrated by Figure 6 in which the error in the measurement of actual viscosity is plotted for the 78 natural 
gases used in Figures 3 -5. 

This invention has been described in this application in considerable detail in order to comply with the Pa- 
tent Statutes and to provide those skflled in the art with the information needed to apply the novel principles 
and to construct and use such specialized components as are required. However, it is to be further understood 
that the invention can be carried out by specifically different equipment and devices and that various modifi- 
cations both as to equipment and procedure details can be accomplished without departing from the scope of 
the invention itself. 
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Claims 

1. Apparatus for determining thermophysical or thermochemical parameters of a fuel gas, the apparatus 
characterised by: 

means (24, 26, 28) to obtain the thermal conductivity, kg, of the fuel of interest at a reference con- 
dition; 

means (24, 26, 28) to obtain the specific heat, Cpg, of a fuel of interest at a reference condition; 
means (24, 26, 28) to obtain dCp/dT at the reference condition of 0,*; 
means (24, 26, 28) to obtain Cpv T1 /Cp, Ts at the reference condition of c^ Ts ; 
means (24, 26, 28) to obtain dk/dT at the reference condition of kg; 

means (36) to determine the desired parameter of the fuel according to a relationship selected from 
Y S = AO + Al*s nl Cps nl *S Pl y S ql T s rl P s 81 + 



+ A 2 k s n2 c p m2 x s P 2 y s « 2 T s r2 P s « + . . . 



— **i A x K s c ps x s Ys T s p s 
where 

Y a represents a thermophysical or thermochemical parameter of interest selected from but not lim- 
ited to any one of higher heating value, H; oxygen demand, D; Wobbe Index, Wo; relative density or spe- 
cific gravity, p ; absolute density p a ; inerts, I; compressibility factor, Z; critical compression ratio, Rc; vis- 
cosity, t] ; 

Ao, A, ... A, are constants or coefficients, 

n i( m,, pi, q,, r b s, are exponents of values from - 20 to 20, including zero, with the number of terms, 
i, ranging from 1 to 15, 

k s and Cpg represent their values at a predetermined reference condition of temperature and pres- 
sure, 

x s represents dk/dT at the reference condition, 

y 9 represents dCp/dT or c^fc^ (T, * T 8 ) at a reference condition, 

T s represents temperature at the reference condition, 

P s represents absolute pressure at the reference condition. 

2. Apparatus according to Claim 1, characterised by means (36) to fit the relation individually to a series of 
limited ranges of values of the parameter of interest which include two or more ranges, for any of the in- 
volved parameters to thereby optimize the relation for each given range. 

3. Apparatus according to Claim 1 or 2, characterised by means (36) to relate the variable Y s to the oxygen 
demand of the fuel D and means to control the firing of one or more burners at a set point based on an 
amount of excess oxygen. 

4. Apparatus according to any preceding Claim characterised in that the parameter of interest is the higher 
heating value. 

5. Apparatus according to any preceding Claim characterised by means (36) to derive the absolute density 
from the values of Cp and from the relation 

Pa = PoVmc/V m 

where 

p 0 is the density of the air at standard conditions (0°C and 1 Atmosphere), 
p A is the density at actual conditions of P and T, absolute specific gravity, 
Vmq is a constant equal to the molar volume at standard temperature air pressure, i.e. » 22-4 l/m 

and 

V M = Cp/Cpv is the actual molar volume at the temperature and pressure of measurement 
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Apparatus according to any preceding Claim characterised in that the parameter of interest is compres- 
sibility, Z. 

Apparatus according to Claim 6 characterised by means (36) to determine Zm B and to correct that value 
by applying the difference 

AZ = Z MB (P,MB.Ti;MB)-2(P r ,T r ) 

where 

Pr = P/Ppo 

T r = T/Tpc 

(Reduced temperature and pressure MB indicates derived from measurement). 

Apparatus for determining thermophysical or thermochemica! characteristics of a fuel gas of interest, the 
apparatus characterised by: 

(a) means (1 6) to conduct at least a partial stream of the fuel gas through a sensor chamber (1 7) having 
one or more sensors in contact with the fuel gas; 

(b) means (24, 26, 28) to generate a first electrical signal indicative of temperature of the fuel gas; 

(c) means (24, 26, 28) to generate a second electrical signal at a sensor in contact with the fuel gas, 
the first electrical signal being representative of a first fuel gas quality, the first fuel gas quality com- 
prising one of the thermal conductivity and specific heat; 

(d) means to conduct the first electrical signal to a computing means (36); 

(e) means (24, 26, 28) to generate a third electrical signal at a sensor in contact with the fuel gas, the 
third electrical signal being representative of a second fuel gas quality, the second fuel gas quality com- 
prising one of the thermal conductivity and specific heat, the second gas quality comprising a different 
one of the fuel gas qualities than the first gas quality; 

(f) means to conduct the second electrical signal to the computing means (36); 

(g) means to repeat steps (a)-(f) at a different fuel temperature; and 

(h) the computing means is to derive a signal for at least one of measurement and regulation using the 
received first and second electrical signals representative of the first fuel gas quality and second fuel 
gas quality, respectively, as a measure for at least one of the desired parameters according to a rela- 
tionship 



+ A 2 k s a2 Cp m2 x s P 2 y a <J2 Ts "p s S2 + 



where 

Y s represents a thermophysical or thermochemica! parameter of interest selected from but not limited to 
any one of higher heating value, H; oxygen demand, D; Wobbe Index, Wo; relative density or specific grav- 
ity, p ; absolute density p a ; inerts, I; compressibility factor, Z; critical compression ratio, viscosity.ti ; 
Ao, A t . . . Aj are constants or coefficients, 

n,, rrii, p,, q lt r,, s, are exponents of values from -20 to 20, including zero, with the number of terms, 
i, ranging from 1 to 15, 

ks and c^, represent their values at a predetermined reference condition of temperature and pres- 
sure, 

x, represents dk/dT at the reference condition, 

y 3 represents dc,/dT or c^/Cp^ (T t * T 8 ) at a reference condition, 

T, represents temperature at the reference condition, 

P, represents absolute pressure at the reference condition. 

Apparatus according to any preceding Claim characterised by means (36) to determine and by the 
expressions: 
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K = ao + a^i - aJte - k^O", - T S )/(T 2 - T,) 
Cp, = b 0 + btCJ* - batCpQ - Cpt1 ) (T, - T 8 )/(T 2 -T t ) 

where 

aj and b ( are constants 

ml and n1 are exponents 

Ti and T 2 are temperatures in a desired range 

T 8 is the reference temperature. 

10. Apparatus according to any preceding Claim, characterised by means (36) to convert each derived signal 
Y into an electrical signal in the form of a fuel gas regulation signal. 

11 Apparatus according to any preceding Claim, characterised by means (36) to automatically periodically 
repeat the determination of Y 8 for the parameter of interest 

12. Apparatus according to any preceding Claim, characterised by means (36) to transmit the derived signals 
for Y s of the parameter of interest to at least one of display means and recording means. 

13. A method of determining thermophysical or thermochemical parameters of a fuel gas, the method char- 
acterised by the steps of: 

obtaining the thermal conductivity, kg, of the fuel of interest at a reference condition; 
obtaining dCp/dt at the reference condition of Cpsi 
obtaining Cp^/Cpy-ft at the reference condition of Cp v7s ; 
obtaining dk/dT at the reference condition of k 8 ; 

determining the desired parameter of the fuel according to a relationship selected from 
Y S = A c + A 1 k s nl c p3 ml Xs^ygql T 3 rl P s sl + 



+ A 2 k 3 *%* 2 x s P 2 y s <I 2 T3 r2 P s 32 + . . . 



- SiAi*s ni c ps ffii *s pi Ys qi T s ri P s si 
where 

Y 8 represents a thermophysical or thermochemicaJ parameter of interest selected from but not lim- 
ited to any one of higher heating value, H; oxygen demand, D; Wobbe Index, Wo; relative density or spe- 
cific gravity, p ; absolute density p a; inerts, I; compressibility factor, Z; critical compression ratio, vis- 
cosity, r\ ; 

Ao, A t ... A, are constants or coefficients, n,, m,, p,, q,, r, s, are exponents of values from - 20 to 
20, including zero, with the number of terms, i. ranging from 1 to 15, 

k, and represent their values at a predetermined reference condition of temperature and pres- 
sure, 

x 3 represents dk/dT at the reference condition, 

y 8 represents dCp/dT or c^n/c^rs (T t * T 8 ) at a reference condition, 

T 8 represents temperature at the reference condition, P 8 represents absolute pressure at the ref- 
erence condition. 

14. A method according to Claim 1 3 characterised in that the relation is fitted individually to a series of limited 
ranges of values of the parameter of interest which include two or more ranges, for any of the involved 
parameters to thereby optimize the relation for each given range. 

15. A method according to Claim 13 or 14 characterised in that the additional steps of: 

relating the variable Y 8 to the oxygen demand of the fuel D; and 

controlling the firing of one or more burners at a set point based on an amount of excess oxygen. 

16. A method according to any of Claims 13 to 1 5 characterised in that the parameter of interest is the higher 
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heating value. 

17. A method according to any of Claims 13 to 16 characterised in that the absolute density is derived from 
the values of C p and Cp V from the relation. 

Pa - PoVmo/v m 

where 

p 0 is the density of the air at standard conditions (0°C and 1 Atmosphere), 
p A is the density at actual conditions of P and T, absolute specific gravity, 
Vmq is a constant equal to the molar volume at standard temperature air pressure, i.e. * 22.4 1/m 

and 

V M = V 0 !" is the actual molar volume at the temperature and pressure of measurement 

18. Amethod according to any of Claims 13 to 17 characterised in that the parameter of interest is compres- 
sibility, Z. 

19. A method according to Claim 18 characterised in that the Z determined is Zm B and wherein this value is 
corrected by the further step of applying the difference 

AZ = Zmb^mbJomb) -ZfPrJJ 

where 

Pr=P/Ppc 
Tr^T/Tpc 

(Reduced temperature and pressure MB indicates derived from measurement). 

20. A method of determining thermophysical or thermochemical characteristics of a fuel gas of interest com- 
prising the steps of: 

(a) conducting at least a partial stream of the fuel gas through a sensor chamber (17) having one or 
more sensors (24, 26, 28) in contact with the fuel gas; 

(b) generating a first electrical signal indicative of temperature of the fuel gas; 

(c) generating a second electrical signal at a sensor in contact with the fuel gas, the first electrical signal 
being representative of a first fuel gas quality, the first fuel gas quality comprising one of the thermal 
conductivity and specific heat; 

(d) conducting the first electrical signal to a computing means (36); 

(e) generating a third electrical signal at a sensor in contact with the fuel gas, the third electrical signal 
being representative of a second fuel gas quality, the second fuel gas quality comprising one of the 
thermal conductivity and specific heat, the second gas quality comprising a different one of the fuel 
gas qualities than the first gas quality; 

(f) conducting the second electrical signal to the computing means (30) ; 

(g) repeat steps (a) - (0 at a different fuel temperature; and 

(h) using the computing means (36) to derive a signal for at least one of measurement and regulation 
using the received first and second electrical signals representative of the first fuel gas quality and sec- 
ond fuel gas quality, respectively, as a measure for at least one of the desired parameters according 
to a relationship 



+ A 2 k s n 2 Cp *2 K 3 P 2 y 3 q 2 T 3 "P s 32 + . . . 



= ZiAi*s ni c p3 mi x s pi y 8 qi T 3 ri P 3 s i 
where 

Y 9 represents a thermophysical or thermochemical parameter of interest selected from but not lim- 
ited to any one of higher heating value, H; oxygen demand, D; Wobbe Index, Wo; relative density or spe- 
cific gravity, p ; absolute density p a ; inerts, I; compressibility factor, Z; critical compression ratio, vis- 
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cosity, ; 

Ao, . . . A| are constants or coefficients, t\ h mi, p u q l( r i( s,, are exponents of values from -20 to 
20, including zero, with the number of terms, i, ranging from 1 to 15, 
x s represents dk/dT at the reference condition, 
5 y 8 represents dCp/dT or o^^fc^ (T, * T 8 ) at a reference condition, 

T 8 represents temperature at the referene condition, 
P 8 represents absolute pressure at the reference condition. 

21. A method according to any of Claims 13 to 20 characterised in that Kg and Cpg are determined by the ex- 
10 pressions: 

ks = ao + a,^ 1 - a^ke - l^fT, - TJ!(T 2 - T, ) 
Cpe = b c + btV - b^Cpc - Cpt1 ) (T, - T 8 )/(T 2 - T, ) 

where 

a, and b { are constants 
15 ml and n1 are exponents 

T i and T 2 are temperatures in a desired range 
T s is the reference temperature. 

22. A method according to any of Claims 13 to 21 characterised in that further comprising the steps of. 

20 converting each derived signal Y into an electrical signal in the form of a fuel gas regulation signal. 

23. A method according to any of Claims 1 3 to 22 characterised by the step of automatically periodically re- 
peating the determination of Y 8 for the parameter of interest. 



24. A method according to any of Claims 13 to 23 characterised by the step of: 
transmitting the < 
and seconding means. 



25 transmitting the derived signals for Y 8 of the parameter of interest to at least one of display means 
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